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Abstract 

The Wheeler-DeWitt equation is derived from the bosonic sector of the heterotic string 
effective action assuming a toroidal compactification. The spatially closed, higher di- 
mensional Friedmann-Robertson- Walker (FRW) cosmology is investigated and a suitable 
change of variables rewrites the equation in a canonical form. Real- and imaginary-phase 
exact solutions are found and a method of successive approximations is employed to find 
more general power series solutions. The quantum cosmology of the Bianchi IX universe 
is also investigated and a class of exact solutions is found. 
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1 Introduction 



Superstring theory remains a promising candidate for a complete theory of quantum gravity 
If string theory is indeed correct, it should describe the complete quantum history of the 
universe. However a consistent second quantization of the superstring is not currently available 
and to make progress one must consider the low-energy limit of the effective field theory. One 
approach adopted by a number of authors has been to investigate the quantum cosmology of the 
superstring by solving the corresponding Wheeler-DeWitt (WDW) equation for the effective 
action [2,3]. 

String theory and the WDW equation may be related at a fundamental level. It has been 
shown that the effective action, S, may be written as an integral over a spatial hypersurface 
that does not depend on cosmic time This is significant because the WDW equation may be 
viewed as a time- independent Schrodinger equation and its solution in the WKB approximation 
is determined by the phase factor exp{iS). Consistency therefore suggests that the action should 
not depend on cosmic time and the above approach appears well motivated. 

Recently the Scherk-Schwarz |Q dimensional reduction technique has been employed to con- 
struct a string effective action in which d internal dimensions are compactified onto a torus [^]. 
In this paper we search for solutions to the WDW equation derived from this action. Broadly 
speaking a singular, oscillating solution may be interpreted at the semiclassical level in terms 
of Lorentzian four-geometries, whereas an exponential solution can represent a quantum worm- 
hole if appropriate boundary conditions are satisfied. Hawking and Page [0 have proposed that 
the wavefunction should be exponentially damped at infinity and obey a suitable regularity 
condition when the spatial metric degenerates (i.e. its radius a vanishes) [§]. Presently the 
precise form of such a condition is not clear, but it appears that the wavefunction should vary 
as \1/ oc a^, where the exponent p is a constant. The sign of p is associated with the choice of 
factor ordering and may be negative 0. 
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The derivation of the string effective action is summarized in Sect. 2 and the correspond- 
ing WDW equation for the spatially closed, [D + l)-dimensional Friedmann-Robertson- Walker 
(FRW) space-time is derived. Singular and non-singular solutions to this equation are found in 
Sect. 3. The dimensionally reduced action is invariant under a global 0{d,d) transformation 
and in principle this symmetry may be employed to generate new solutions. Such a procedure 



is illustrated within the context of the duality symmetry of string theory |]10l . The WDW equa- 
tion is rewritten in a simplified canonical form after a suitable redefinition of the independent 
variables and this allows a method of successive approximations to be developed in Sect. 4. 
Additional power series solutions are found. Some of these series can be written in closed form 
and are exponentially damped at infinity. The more general Bianchi IX (mixmaster) cosmology 
is investigated in Sect. 5 and an exact solution is found that is sharply peaked around the FRW 
universe at large three-geometries. 

Units are chosen such that Ti = c = IGttG = 1. 

2 The Effective Action and the Wheeler-DeWitt Equa- 
tion 

We first summarize the derivation of the string effective action presented in At the tree 
level, the Euclidean action of the gauge singlet, bosonic sector of the D-dimensional heterotic 
string isQ 



S-g= d^x\lge 

J M 



(2.1) 



where (f) is the dilaton field and H is the totally antisymmetric 3-index field [|TT|. In the technique 

of dimensional reduction the universe is viewed as the product space Ai = x JC, where the 

[D + l)-dimensional space-time J^{x^) has metric g^y{xP) and the ci-dimensional internal space 
hat "denotes quantities in the ZJ-dimensional spacetime Ai, where D = 1 + D + d. 
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/C(y") must be Ricci flat if tlie matter fields are independent of its coordinates y". This is the 
case for the Calabi-Yau spaces often considered in string theory, but for the purposes of the 
present work it is sufficient to assume /C is a torus, i.e. )C = x x . . . x . The complete 
metric is then 

where Gajs is the metric on /C. The effective action in {D + l)-dimensions is 

S = J d^'+'x^e-^ [-i? - (V0)2 + ^H^^pH^^'P - ^Tr (v^ikf- V^M 



where 



(2.3) 



= + B^,AW^ (2.4) 
= 0--lndetG (2.5) 

Hf,up = V^,B^p - ^A'^r]ijJ^lp + (cyc.perms.) (2.6) 
= - V.4 (2.7) 

and the 2d x 2d matrices M, t] and are defined by 

respectively. It follows from the definition of M that M^rjM = rj and this implies M G 0{d,d). 
Hence, since g^u and (p are also invariant under the action of this group, the dimensionally 
reduced action is symmetric under a global 0{d,d) transformation. 

We follow the approach of by choosing = constant with H^^p = and = 0. We 
assume a line element for J' of the form 

ds"^ = n^{t)dt^ + a^{t)dnl, (2.9) 



where n{t) defines the lapse function, dQ% is the line element of the unit D-sphere and a{t) is 
the scale factor. This assumption of spherical symmetry implies that all variables are functions 
of t only. Thus the d x d matrix G + B may be written as 



G' + 5 = diag(Ei,...,Ed/2), S,= (^ 



(2.10) 



where G and B are space-time dependent and we assume d is even. The action ( p.3| ) then 
simplifies to 

/ . d/2 \ 

S = j rf^+^xV^(-i? + -E[(V^.)' + e-'''^(V^.)'] 1 • (2-11) 

One special example of an 0{d,d) invariance is the duality symmetry associated with an 
interchange of M and [jlO[. The duality transformed fields are 

(e'^. +e-^^a|)~'e-^^(T,- (2.12) 



a,- = - 



and a new class of solution parametrized by {V'jj'^j} i^^iy be generated from a solution given 
in terms of {^/'j, a^). 



The WDW equation for the theory ( 2.11 ) may now be derived using the standard techniques 



12 . The classical Hamiltonian constraint is 



d/2 

a^n oc -a^n^ + 2D{D - 1) ^ [h^^, + e^'^^H^J - ^^^^(d-i) ^ ^2.13) 

where 7i is the Hamiltonian density, IIx. = dC/dzi are the momenta conjugate to Zi and a dot 
denotes differentiation with respect to t. The constant a = 2D{D — l)Qu where i^/) is the 
volume of the unit -D-sphere. The system is quantized by identifying the conjugate momenta 
with the operators 
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where the arbitrary constants {p,qj,rj} account for ambiguities in the operator ordering. The 
WDW equation follows by viewing 7i as an operator acting on the wavefunction \E'[a, ipj, aj]: 



o o 

oa oa 



d/2 I 

i=i V 



(T ■ 



-a 



We proceed in the following sections to solve this equation. 



3 Exact Solutions 



(2.15) 



Exact solutions to Eq. ( p.l5| ) may be found for qj = 0. We first search for separable solutions 
of the form 



d/2 

The WDW equation decouples to {d + 2)/2 differential equations: 



(3.1) 



d „ d 



da da 



,=1 \di^] 



$ = 



(3.2) 



a, 



aY^+uj^Q = Q, 1 = 1,..., d/2, 



d(Ti dai 

where ui are arbitrary separation constants. Eq. ( p. 3D admits the exact solutions 



Ci = constant, ui = 0, 



(3.3) 



(3.4) 
(3.5) 



where Jp is a Bessel function of order p. Classically Eq. (|3.5|) corresponds to d/ = 0, which is 
a particular solution to the Einstein field equations. In this case each of the cx/ may be set to 
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zero without loss of generality by means of a linear translation. Thus we shall first solve Eq. 
( p. 21 ) for uji = 0. In this paper we are concerned only with the functional forms of the solutions 
to Eq. ( [^.15| ) and we therefore ignore all constants of proportionality in the wavefunction. 



3.1 Case A: = 



It proves convenient to make a change of variables to {u, v,Si, . . . , S(rf_2)/2} defined in such a 
way that <l> is independent of the Sj. This necessarily results in some loss of generality but 
allows analytical solutions to be found. It is shown in the appendix that with a factor ordering 
p = I, Eq. ( |3.2|) reduces to the simplified canonical form 



d d 

du dv 



a 
T 



\uv 



:D-2 



$ = 0, 



where 



u = a exp 



2DiD - 1) ,=1 



d/2 

^3 COS Oj 



(3.6) 



(3.7) 



V = a exp 



d/2 



2D{D - 1) ,=1 



(3. 



and the set of constants {Qj\ are solutions to the constraint equation 



d/2 



(3.9) 



The {u,v} variables are useful because they eliminate the direct dependence of the WDW 
equation on the number of axion fields present. We shall refer to them as the canonical coordi- 
nates. These variables are the null coordinates over minisuperspace when only one axion field 
is present [|13[ . The form of Eq. (p.6|) remains invariant if one adds or removes extra massless 



scalar fields and the alterations are accounted for by changing the number of fields included 
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in the summations of Eqs. (|3] 
interchange of these variables. 

The variables u and v may be redefined as 



]9D. We note also that Eq. ( p.6|) is invariant under an 



u 



(— 

V 2a 



1/{D-1) 



In this case Eq. ( |3.6| ) transforms to 
92 



2 I 2 
y + X 



$ = 0, 



dx"^ 

which is the equation for two harmonic oscillators with equal and opposite energy 
A singular solution to Eq. ( |3.2| ) is 



(3.10) 



(3.11) 



$^ = K±i, (aa^-^/{D - 1)) exp 



d/2 

Z7 ^ ipj COS Qj 



where 



d/2 

cos^ 9 J = 1, 

i^' is a modified Bessel function, 7 is an arbitrary constant and 



7_ D-l 

7 = V^' 



On the other hand, non-singular solutions to Eqs. ( |3.6| ) and ( |3.r 



are 



$c = exp 



2(L) 



3n (^-'''"-' 



D-l 



and 



•^-n = ^if„(a;)//„(i/) exp [- (x^ + y^) /2 



(3.12) 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



respectively, where c is an arbitrary integration constant and iJ„ are Hermite polynomials 
Solution ( p.l5|) may be written directly in terms of the scale factor and axion fields: 



a 



A = exp 



D - 1 



cosh 



'd/2 



- I EV^iCos^j I + A 



(3.17) 



where A = Inc. Solutions equivalent to ( p.l6|) and (|3.17|) have been found previously in a 



different context by Zhuk ||T^ . They do not oscillate as the spatial metric degenerates and they 



are exponentially damped at infinity. Consequently they may be viewed as quantum wormholes 
with asymptotic topology 3? x S^. 

The ground state wavefunction of Eq. ( ^.151 ) corresponds to c = 1 (A = 0) and excited states, 
corresponding to c 7^ 1, are generated by linear translations on the ipj fields [If. Thus there 
exists a continuous spectrum of these excited wormhole states. On the other hand Eq. ( 3.16| ) 



represents a descrete spectrum of excited states 0, but the two spectra have the equivalent 
ground state. The duality transformation ( p. 12 ) maps ipj to —ipj when the aj fields vanish 



and this is equivalent to an interchange of the canonical coordinates. Hence it follows from 
Eq. (|3.15|) that the ground state may be interpreted as the self-dual solution to the WDW 



equation. In this sense the ground state is associated with the solution of maximum symmetry. 
In figure 1 this ground state and an excited state of the continuous wormhole spectrum ( |3.15| ) 
are plotted to illustrate this feature. 

Figure 1 

Additional wormhole solutions may be found by factoring out the continuous spectrum given 
by Eq. (|3.15| ). If we write 

v) = A{u, v)<^c{u, v), (3.18) 

it follows after substitution into Eq. ( p.6|) that Eq. (|3.18|) is a solution to the WDW equation 
if A satisfies 

d'^A acv^-^ OA au^-^ - p (3 I9) 

dudv 2 du 2c dv 

When c = 1 a comparison with Eq. (|3.16|) implies that A may be written as a product of two 
Hermite polynomials. For c 7^ 1 one solution to Eq. ( p.l9| ) is 

A = ct;^-^ - c-^M^-^ (3.20) 
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2a^-^ sinh 



7 



'd/2 



(3.21) 



3.2 Case B: ui ^ 



The question now arises as to whether singular and non-singular solutions to Eq. ( |3.2| ) can be 
found when ui ^ 0. If we assume the separable ansatz 



d/2 

$[a,^,] = A(a)ni?,(^.) 
for arbitrary functions A{a) and Bj{ipj), we find 

^ / ^ _ 2 _ ^,2^2(^-1)) A = 

da an. V ^ ' / 



and 



1 (PBj 2 



2 = _^^^ = _^^cos^e„ 



(3.22) 



(3.23) 



(3.24) 



where Qj again satisfy the integrability condition ( |3.13| ) and z is an arbitrary constant. 

Eq. ( ^.241 ) reduces to a Bessel equation after the change of variables ipj = ^^^j- Hence one 
solution to Eq. ( 3.2|) is 



d/2 



= a^'-^y'K^s (aa^-V(I) - 1)) [] {^^e^ 

i=i 



(3.25) 



where 



(D-l) 



1 — p 



- 2D{D - l)z^ 



(3.26) 



and the expression for the full wavefunction \1/ is 



d/2 



vl>[a, V',, a,] = a^'~^^/^K, (aa^- V(D - 1)) 11 (^^e^^O ^^^'^i) ^r^'^^' • (3-27) 
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For simplicity let us consider the case p = 1 and rj = 0. This leaves the WDW equation 
invariant under a change of minisuperspace coordinates |13[ and the Cj functions represent 
plane wave solutions. The modified Bessel function Kq{x) is exponentially damped at infinity, 
but diverges as Kg oc x~'^ in the limit x — 0. It follows from Eq. ( |3.25| ) that 



$z oc exp 



—tz 



d/2 



2D{D - 1) In a + ^ cos 9^- 



(3.28) 



as {a, e^^} and the wavefunction oscillates an infinite number of times. This corresponds to 
an initial singularity, but may be removed by integrating over the separation constant z. This 



is analogous to taking the Fourier transform p|. The Riemann-Lebesgue lemma states that 



lim 



dxe'P'' = 



(3.29) 



and it follows that the transformed wavefunction 



+ 00 



dz^:. 



(3.30) 



is damped both at the origin and at infinity. It is therefore Euclidean for all values of a. 

Finally new solutions may now be generated from Eq. ( p.27| ) in terms of the duality trans- 
formed fields defined in Eq. ( |2.12| ). After a simple rearrangement we find 



-1 



(3.31) 



and 



(1 



(3.32) 



Substitution of these expressions into Eq. ( p. 27] ) generates a more complicated class of singular 
solution [a, ipj , 5"j] . 
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4 Power Series Solutions via a Method of Successive Ap- 
proximations 



Additional solutions to Eq. ( p.6| ) may be derived by employing a method of successive approx- 
imations. We assume the wavefunction may be expanded as the infinite sum 

oo 

= J2 ^m{u,v). (4.1) 
m=0 

This is a solution to Eq. ( |3.6| ) if the functions $o and {m > 1) satisfy 

(4,2) 



dudv 
and 

4 



respectively. Eq. ( [4.2|) is the canonical form of the one-dimensional wave equation and admits 



the general solution 

= f{u) + g{v), (4.4) 

where f{u) and g{v) are arbitrary twice continuously differentiable functions. 

An iteration procedure may now be established by substituting the solution for $0 into Eq. 



( [4. 3D , solving for $1, and then repeating the process. The symmetry of Eqs. (|4.3| ) and 
allows this to be done analytically by expressing the solutions in terms of quadratures with 
respect to the canonical coordinates. For simplicity we assume for the moment that g{v) = 0. 
It follows that for m > 1 niay be written in terms of quadratures with respect to u: 

*m=(^-j rn\{D-l)"' J dUm^iuZ~l ■ ■ ■ j duiu^-"^ f {ui). (4.5) 

By symmetry however, if $0 = g{v), the form of is identical to Eq. ([4.5|) but with u replaced 
by V and f{ui) replaced by g{vi). A more general expression for the wavefunction is therefore 
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Solution 



I 

II 
III 



D-l 



l-D 



Wavef unction 



D-l 



U 



it; I sinh 



D-l 



cosh 



(md) 



{D-l)/2 



Table 1: Closed solutions to the canonical form of the WDW equation when $o — v^ztu^. /g is the modified 
Bessel function of the first kind of order zero. All solutions diverge at infinity, but linear combinations may be 
found which are exponentially damped at infinity. It is straightforward to verify by differentiation that these 



are solutions to Eq. (3.6) 



u 



m(D-l) / 



dUm-lU. 



D-2 
m—l 



U2 



dUiU^ f{Ui) 



(4.7) 



Gmiv) 



1 



V 



m{D-l) 



dv V^~^ 



dvm-iv'^Ji ... I dviv^ "^givij 



The integrals in Eqs. ( |4.7| ) and ([4.8|) may be evaluated analytically when 



= v'^ ±. u^, b = constant 



and it is straightforward to show that 



F -\- G 



[D-1+ b][2{D - 1) + 6] . . . [m{D -!) + &]' 
It follows that the wavefunction reduces to the series solution 

1 



$ = <l>n 



1+E 



oo / X 9jTi / \m(D—l) 



\\2J m!(D - 1)™ [D - 1 + 6] . . . [m{D - 1) + 6] 



(4.9) 



(4.10) 



(4.11) 



Closed expressions for $ exist when 6 = and h = ±{D — l)/2 and these solutions are 
shown in Table 1. They are Euclidean for all values of the scale factor in the sense that they do 
not oscillate. However they diverge exponentially fast at infinity and therefore do not represent 
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quantum wormholes in the Hawking-Page sense. Solution I is identified as the Hartle-Hawking 



ground state [0. On the other hand a hnear combination of solutions II and III is 



$111 - $T 



D- 1 



1/2 



Kl/2 



a 



,D-1 



cosh 



I i^i cos 9i 

8D J f^J^ ^ 



(4.12) 



\tt{D-1)J "''^\D-1 
This solution has the correct asymptotic behaviour at infinity to be interpreted as a quantum 
wormhole. Although it diverges as a^^"^^/^ when the spatial metric degenerates, this may be 
acceptable behaviour for a quantum wormhole 0. Moreover, Eq. ( p.6| ) was derived from the 
field theoretic limit of the string effective action, and one might expect the analysis of Sect. 2 
to be invalid for scales below the Planck length. One could therefore argue that this divergence 
will not arise in a more accurate analysis. 



5 An Exact Solution in the Bianchi IX Cosmology 



The dimensionally reduced action ( p.ll| ) is essentially {D + l)-dimensional Einstein gravity 
minimally coupled to a set of massless scalar fields with non-standard kinetic terms. It is 
therefore possible when D = 3 to extend the results from the previous sections to the class of 
anisotropic and homogeneous Bianchi cosmologies. The line element for the class A Bianchi 
spaces is 



ds^ = -dr + e 



2a(t) 



(5.1) 



where e is the determinant of the metric on the three-surface, (3ij is a tracefree, symmetric 
matrix and are the one-forms determining the isometry of the three-surface. Without loss of 
generality [3ij may be diagonalized: 



A,(t) = diag [(3+{t) + y3/5_(t),/?+(t) - v^/5_(t), -2/5+(t) 



(5.2) 
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Following an identical procedure to that presented at the end of Sect. 2, the WDW equation 
derived from the action ( |2.11| ) is found to be 



da da 



dpl 80: 



\ + t/(a,/?±)-12Ehr72 



^ = 0, 



(5.3) 



where the potential ?7(a, depends on the spatial curvature in the model. For simplicity 
we choose a factor ordering qj = r j = for the scalar field momentum operators and rescale 
such that a = 1. The WDW equation ( p. 15 ) for the positively curved, isotropic FRW model 
is recovered by setting /5+ = /3_ = in the potential and removing the d'^'^/df3'^ terms in Eq. 
( ^.3|) . When /3_ = the Bianchi IX universe simplifies to the Taub cosmology ||T8| . 
To proceed we assume the separable ansatz: 



[a, p+, /?_, V-i, a,] = $ [S, p+, S [ipj, aj] . 



The WDW equation separates into two differential equations: 



d 



(9^ ^ ^ ^ 



$ = 



(5.4) 



(5.5) 



d/2 



,=1 V^^J 



12 E 



+ e 



'da]. 



— m 



5 = 0, 



where m is an arbitrary separation constant. 

The curvature potential for the type IX may be written in the form 



3f/ = e 



4a 



e-8^+ - 4e-2^+ cosh 273/3. + 2e^^+ (cosh4y3/5_ - 1 



(5.6) 



(5.7) 



It is well known that Eq. ( ^.7|) exhibits a triangular symmetry in the (/?+,/?_) plane, but 
recently it was shown pO|] that there exists an additional symmetry 

f dxY ( dx 



dx 
da 



dp, 



dp- 



+ U = 0, 



14 



where 

X = -e^^Tr f e^^) = -e^'^ fe-^^+ + 26^^+ cosh 2 Vs/?.! . (5.9) 

g V / 5 L J 

This additional symmetry may be employed to find new exact solutions. If we substitute 
the secondary ansatz 

$ = Vr(S)e-^ (5.10) 



into Eq. (|5.5|) , and employ the identities 

d^X d^X d^X . dx 



aa^ QUI 9/3! + "'^ = at - = »• 

it follows that Eq. ( ^.5|) simplifies to 

d'^W dW dW 

-—r- + 12xW^ + - 2pxW - - m^W = 0. (5.12) 

For specific choices of the factor ordering parameter p, Eq. ( |5.12| ) admits the exponential 
solution W = exp(— ca), where c is a decay constant. Eq. ( ^.12| ) reduces to the constraint 
equation 

-pc-m'^ + 2{p- 2c- 6)x = 0, (5.13) 
which is solved for 



p = 2(c + 3), c= -3±i/(3 + m)(3-m). (5.14) 

Requiring the factor ordering to be real implies |m| < 3. This solution is a generalization of 
the vacuum solution corresponding to m = found in . We find a new vacuum solution for 
c = p = —6. It is worth remarking that the factor in Eq. ( |5.10| ) appears to be independent 
of the choice of factor ordering. 
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Finally we may combine these results with the solution to Eq. ( ^.6| ). This equation is solved 
by separating it into Eqs. ( pl3| ) and ( |3.24| ) as in Sect. 3. Hence, one solution to Eq. ( |5.3| ) for 
the Bianchi IX cosmology is 

^ = e--->^ n {^^ cose, {co,e^^) e^^"^'^^ } , (5.15) 

where Uj are arbitrary constants that may be imaginary and the 6j must satisfy J^fJi cosO'j = 1. 
A new class of solution may now be generated from Eq. (|5.15 ) as a consquence of the 



duality symmetry of the action ( p. 11 ), as summarized in Eqs. ( |3.31 ) and ( 3.32 ). However, the 



real-phase part of the solution, $, remains invariant after such a transformation. 



When c = it was noticed in |£T| that as d increases the wavefunction $ becomes peaked 
in the (/3+,/3_) plane around the isotropic FRW solution /5+ = /5_ = 0. Our solutions also 
exhibit this feature for values of c given by Eq. ( ^.13 ). In principle, a solution of the form 



$ = e~'^"'~^ exists for any WDW equation that contains a separable component of the form 
( p.5| ). For example such a separation is possible with any matter sector that behaves as a stiff 
perfect fluid at the classical level. Dynamically this is equivalent to a massless scalar field and 
at the quantum-mechanical level the matter is an eigenstate of some operator with eigenvalue 
m. The constant 27r^m represents a conserved scalar flux. 



6 Discussion 



In this paper the WDW equation was derived from the dimensionally reduced string effective 
action assuming a toroidal compactification. The {D + l)-dimensional FRW and the four- 
dimensional Bianchi IX (mixmaster) cosmologies were investigated. In the former example 
singular solutions were found and these may be interpreted as spatially closed Lorentzian uni- 
verses. Euclidean solutions were also found which may be viewed as quantum wormholes if 
appropriate boundary conditions are satisfied. The correct quantum theory of gravity should 
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explain, either from symmetry or probabilistic considerations, why compactification to a four- 
dimensional space-time is observed. For this reason the dimensionality of the external space 
was left unspecified in the analysis. 

In the Bianchi type IX case exact solutions were obtained for specific choices of factor order- 
ing. The wavefunction is strongly peaked around the FRW universe at large three-geometries. 



Recently a symmetry of the form (|5.8|) has been found for the Bianchi type II metric and 
it appears that this symmetry exists for all class A spaces. In principle solutions of the form 
Eq. (|5.15|) can therefore be found for these spaces. 

In both the FRW and Bianchi IX examples the 0{d,d) symmetry of the string effective 
action may be employed to generate new classes of solutions which may otherwise be difficult 
to derive using standard techniques. This was illustrated using the duality symmetry of the 
action, which is a subgroup of 0{d,d). The action ( 2.11| ) is also symmetric under the group 
SL(2, R) and new solutions may also be generated from this symmetry [Q. 

A method of successive approximations was employed for the FRW cosmology to derive 
further solutions. This method should be useful for finding perturbative solutions when uv = a? 
is small (but not so small that the analysis of Sect. 2 breaks down). In this region one may 
view the power series as a truncated Taylor series up to some appropriate order in the scale 
factor. 

Finally, in the regime where the curvature term of Eq. (|3.6|) becomes negligible, the Wheeler- 
DeWitt equation reduces to the one-dimensional wave equation: 

_^ ^^^^ 

dudv 

with general solution $ = f{u) + g{v). Eq. ( |6.1| ) may also be treated as the exact WDW 
equation for a spatially fiat Friedmann universe with nontrivial spatial topology such as the 
3-torus = X X ||2^ . This illustrates a further advantage of employing the canonical 
coordinates. In this case one may easily find forms of f{u) and g{v) which satisfy the Hawking- 
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Page boundary conditions and this suggests that a positive spatial curvature is not a necessary- 
condition for the existence of quantum wormholes ||24|| . 

Acknowledgments The author is supported by the Science and Engineering Research 
Council (SERC) U K and is supported at Fermilab by the DOE and NASA under Grant 
NAGW-2381. 



18 



References 

1. Green M B, Schwarz J H and Witten E 1988 Superstring Theory (Cambridge: Cambridge 
University Press) 

2. DeWitt B S 1967 Phys. Rev. 160 1113 

Wheeler J A 1968 Battelle Rencontres (New York: Benjamin) 

3. Enqvist K, Mohanty S and Nanopoulos D V 1987 Phys. Lett. 192B 327 
Enqvist K, Mohanty S and Nanopoulos D V 1989 Int. J. Mod. Phys. A 4 873 
Pollock M D 1989 Nucl. Phys. B 315 528 

Lyons A and Hawking S W 1991 Phys. Rev. D 44 3802 
Pollock M D 1992 Int. J. Mod. Phys. A 7 4149 
Wang J 1992 Phys. Rev. D 45 412 

4. Giddings S B and Strominger A 1988 Nucl. Phys. B 306 890 
Pollock M D 1991 Nucl. Phys. B 355 133 

5. Scherk J and Schwarz J H 1979 Nucl. Phys. B 153 61 

6. Khastgir S P and Maharana J 1993 Phys. Lett. 301B 191 
Maharana J and Schwarz J H 1993 Nucl. Phys. B 390 3 

7. Hawking S W 1990 Mod. Phys. Lett. A 5 453 
Hawking S W and Page D N 1990 Phys. Rev. D 42 2655 

8. Garay L J 1991 Phys. Rev. D 44 1059 

19 



Kiefer C 1988 Phys. Rev. D 38 1761 
Campbell L M and Garay L J 1991 Phys. Lett. 254B 49 
Wada S 1992 Mod. Phys. Lett. A 7 371 
Gonzalez-Diaz P F 1990 Mod. Phys. Lett. A 5 1307 
9. Alty L J, D'Eath P D and Dowker H F 1992 Phys. Rev. D 46 4402 

10. Kikkawa K and Yamasaki M 1984 Phys. Lett. 149B 357 
Sakai N and Senda I 1986 Prog. Theor. Phys. 75 692 
Pollock M D 1989 Phys. Lett. 227B 221 

Veneziano G 1991 Phys. Lett. 265B 287 

Meissner K A and Veneziano G 1991 Mod. Phys. Lett. A 6 3397 

11. Fradkin E S and Tseytlin A A 1985 Phys. Lett. 158B 316; Nucl. Phys. B 261 1 
Callan C G, Friedan D, Martinec E J and Perry M J 1985 Nucl. Phys. B 262 593 

12. Kolb E W and Turner M S 1990 The Early Universe (New York: Addison- Wesley) 

13. Hawking S W and Page D N 1986 Nucl. Phys. B 264 185 

14. Zhuk A 1992 Phys. Rev. D 45 1192 

15. Bender C M and Orszag S A 1984 Advanced Mathematical Methods for Scientists and 
Engineers, (Singapore: McGraw-Hill) 

16. Hartle J B and Hawking S W 1983 Phys. Rev. D 28 2960 

Zhuk A 1992 Class. Quantum Grav. 9 2029 

20 



17. Ellis G F R and MacCallum M A H 1969 Commun. Math. Phys. 12 108 

18. Taub A 1951 Ann. Math. 53 472 

19. Misner C, Thorne K and Wheeler J A 1973 Gravitation, (Freeman: San Francisco) 

20. Graham R 1991 Phys. Rev. Lett. 67 1381 

21. Moncrief V and Ryan Jr. M P 1991 Phys. Rev. D 44 2375 

22. Obregon O, Socorro J and Benitez J 1993 Phys. Rev. D 47 4471 

23. Grishchuk L P and Sidorov Yu V 1987 Proceedings on the Fourth Seminar on Quantum 
Gravity, (Singapore: World Scientific) 

24. Coule D H and Solomons D 1992 University of Cape Town Preprint 



21 



Appendix 

The purpose of this appendix is to derive the expressions for the canonical coordinates when 



ujj = 0. We make an arbitrary change of variables in Eq. (|3.2| ) and search for solutions that 
depend on only two of these new variables u and v. We choose p = 1 and define a = Ina for 
convenience. The differential momentum operators of Eq. ( p.2|) take the form 

+ hr- + + ittt, (A.i) 



dw"^ dw"^ du \ dw I dv? dw dw dudv dw"^ dv 



where the variable w is to be identified with Ina and ipj. 

One may arrange for the majority of these terms to cancel after substitution into Eq. 
Only the contributions containing d'^^/dvdv remain if u is a solution of the (i/2-dimensional 
wave equation: 

subject to the integrability condition 



Equivalent equations also apply for v. The first integral of Eq. ( A.2|) is 



dv , cos Oj dv / . . N 

- ^ (A.4) 



di^j J2D{D - 1) da ' 



where the constants of integration 9j satisfy the constraint equation Yfj=i cos^ 6j = 1. It follows 
that 



V 



d/2 

f I lna± [2D(D- l)]"^/^^7/^j Cos^j- | , (A.5) 
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where / is some arbitrary function. To ensure that the d'^'^ /dudv terms do not cancel in the 
WDW equation, the argument in v for the scalar fields must take the opposite sign to that in 
M, i.e. 



d/2 

^ I Ina^ [2D(D - 1)]"^/^ ^ t/-^- cos^y | , (A.6) 



where g is a, second arbitrary function. 
If we assume the separable ansatz 

d/2 d/2 

u = u{d) Yl Uj{ipj), V = v{a) Y[ Vjiijj), (A.7) 
j=i j=i 

it is easy to verify by differentiation that one solution to these equations is given by Eqs. (|3.7|) , 
and (Km. 
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Figure Caption 



Figure 1: In Figure (la) a schematic plot of tlie self-dual ground state of both the contin- 
uous and descrete wormhole spectra ( ^.151 ) and ( |3.16D is shown as a function of the canonical 
coordinates. The symmetric nature of the solution is clearly illustrated. In Figure (lb) an 
excited state of ( p.l5|) is plotted for c = 1/2. As noted in the text the ground state may be 
viewed as the solution of maximum symmetry. 
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This figure "figl-l.png" is available in "png" format from: 
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